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Abstract 

Terms contributing to the hyperfine structure of the muonium ground state 
at the level of few tenths of kHz have been evaluated. The a 2 (Za) radia- 
tive correction has been calculated numerically to the precision of 0.02 kHz. 
Leading ln(Za) terms of order a 4 ~ n (Za) n ,n = 1,2,3, and some relativistic 
corrections have been evaluated analytically. The theoretical uncertainty is 
now reduced to 0.17 kHz. At present, however, it is not possible to test QED 

to this precision because of the 1.34 kHz uncertainty due to the muon mass. 
PACS numbers: 36.10.Dr, 12.20.Ds, 31.30.Jv, 06.20.Jr 
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The hyperfine splitting of the muonium ground state is one of very precisely measured 
quantities 

Az/(exp) = 4 463 302.88 (16) kHz (0.036 ppm). (1) 

Currently new experiment is in progress to improve the measurement of Ai/(exp) and muon 
mass by a factor of five or more 0. This is very important for testing the validity of quantum 
electrodynamics (QED) since Au can be calculated very precisely in QED, being relatively 
free from the effect of hadronic interaction. The precision of such a test is limited at present 
by the uncertainty in theoretical calculation, which may exceed 1 kHz. This paper reports 
our result in which we have reduced this uncertainty by nearly an order of magnitude. 
As is well known, the bulk of the hyperfine splitting is given by the Fermi formula 
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where Z is the charge of the muon in units of the electron charge, is the Rydberg constant 
for infinite nuclear mass, and m e and are the electron and muon masses, respectively. 
Of course Z = 1 for the muon, but it is kept in the formula in order to distinguish the 
contribution of binding effect (Zee) from that of radiative correction (a). 

Many correction terms of both a and Za type have been calculated over 40 years. It 
is customary to classify them into three types: radiative non-recoil correction, pure recoil 
correction, and radiative-recoil correction. In addition there is a small weak interaction 
contribution. Thus one may write 

Az/(theory) = Az/(rad) + Az/(recoil) + Az/(rad- recoil) + Az/(weak). (3) 

Conventionally, the effect of hadronic vacuum polarization is included in Ai/(rad-recoil). 

Purely radiative terms of orders a(Za) and a(Za) 2 have been known for some time 0: 
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Here a e and a M are the anomalous magnetic moments of the electron and muon, respectively. 
The appearance of the factor (1 + a M ) in (Q) is in accord with our definition of Ep in ([|). 
The known recoil corrections add up to [Rf 
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where 7 = Zam r ,m r = m e m M / '(m e + m^). The radiative-recoil contributions, which arise 
from both lepton lines and vacuum polarizations, are given by 
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The a{Za) term is known exactly [Q,|J except for the hadronic vacuum polarization contri- 
bution (the (2.15±0.14) term) Q. The In 3 and In 2 parts of the a 2 (Za) term were evaluated 
by Eides et al. 0. Finally there is a small contribution due to the Z° exchange. Our 
re-evaluation of the standard-model estimate M gives 



Az/(weak) ~ 0.065 kHz. 
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As is clear from these results one must know the a 2 (Za) radiative correction in order 
to improve the theoretical prediction further. Fig. 1 shows typical diagrams contributing 
to this order. Recently, terms represented by the diagrams (a) - (e) of Fig. 1 have been 
evaluated by Eides et al. ||. Their results are as follows: 
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= -0.369 kHz, (10) 



a 2 (Za) 

Ai/(Fig.l(d)) = -0.310 742 •• • 1 ' E F 

= -0.171 kHz, (11) 



where z = ln((l+v5)/2). The results (||), © and (0) are analytic, while (|TT|) was evaluated 
numerically after reducing the integral to one dimension. We confirmed these results by an 
independent numerical calculation. However, our purely numerical evaluation of Fig. 1(e): 

Az/(Fig.l(e)) = -0.472 48 (9) a ^——^- Ep 

7T 

= -0.261 kHz (12) 



disagrees with the semi-analytic result of Ref. ||. Recently Eides |IIJ found an error in the 
Table after Eq. (23) of Ref. ||. Their corrected value is in good agreement with (|I2"D . 

Fig. 2 shows the complete set of Feynman diagrams of type (f ) of Fig. 1, which has not 
yet been evaluated. The primary purpose of this paper is to report a preliminary result of 
our calculation for all diagrams of Fig. 2 carried out in the Feynman gauge: 

a 2 (Za) 

Ai/(Fig.l(f)) = (-0.63 ± 0.04) 1 > E F 

= -0.347 (0.022) kHz, (13) 

where the error is mainly due to the uncertainty in extrapolating the integral to zero infrared 
cutoff. Details of calculation will be reported elsewhere. The complete a 2 (Za) correction is 
the sum of ® — (0), (03), and ([13]): 

Ai/(Fig.l) = 0.449 (0.022) kHz. (14) 



Recently we have received two preprints from Eides and his collaborators [JTTJ] , which 
report the result of their calculation, carried out in the Fried- Yennie gauge, for part of the 
diagrams of Fig. 2. We were able to compare our result for the sum of diagrams H17, H18 



and H19 with theirs since this sum is invariant under the covariant gauge transformation. 
They are in perfect agreement although individually they have quite different values due 
to different gauges. To compare other diagrams, we have to wait for completion of their 
calculation. 

The remaining theoretical uncertainty in Av comes mainly from terms of orders a 4 
and a 3 (m e /m^). Although these diagrams are of higher order than ([14]), they may have 
numerically comparable magnitudes due to the appearance of ln(Za) and/or ln(m M /m. e ) 
factors. Some of these contributions are known: [|12] 
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= 0.287 kHz, (15) 

where the first three terms come from the magnetic form factor correction to the 5-function 
potential Vp whose expectation value is Ep, the reduced mass correction to the ln 2 (Za;) 
terms of (f|), and the In k part of the Salpeter term of the Lamb shift | 13| , respectively. The 
last term is a higher order Breit correction. 

Here we report additional terms evaluated in the NRQED perturbation theory ||14j| : 

Au(a) = < V 50 GVp >+••■, 
Av(P) = <V 2 ^ oop GV F >+■••, 
Ai/( 7 ) = < V s GVp >+•••, 
Au(5) = < V Ms G(KotD) >+•••, 

Au(e) = a e < VioGVp >+••■, (16) 

where < • • • > means the difference of the triplet and singlet expectation values with respect 
to the non-relativistic wave function of the muonium ground state. G is the Green's function 
of the non-relativistic electron in the muon Coulomb potential, and D and K are the Darwin 
and /c 4 -kinetic energy term. V50 and V2-100P are 5-function potentials whose expectation 



values are the A50 term of the Lamb-shift energy with one virtual photon spanning over 
any number of Coulomb photons [15j and the term arising from two spanning photons 
as well as the vacuum-polarization effect [IE], respectively. Vs is the 5-function potential 
corresponding to the part of the Salpeter term not included in flT5|). Vmb is the effective 
hyperfine interaction potential corresponding to the In {m^/m e ) term of ([S]). The last term 
of ( |l6l) comes from the A^q part of the Lamb shift and the electron anomaly correction 
to V F . 

All these contributions have been evaluated analytically. Their numerical values are 



Au(a 
Av{(3 

Au(5 
Av{e 



-0.381 kHz, 
-0.007 kHz, 
-0.189 kHz, 
-0.210 kHz, 
0.004 kHz. 



(17) 



The term Au(a) is of order a(Za) 3 ln(Za) and was obtained by Lepage |T7|. The term 
Av(f3) is of order a 2 (Za) 2 ln(Za). Unfortunately, this evaluation is incomplete since V2_i 00 p 
is not yet fully known. However, the contribution of remaining terms will not be much larger 
than the above result. The term Az/(7) is of order {m e /m^){Za) 3 \\i{Za). Au(S) and Az/(e) 
are proportional to {Za) 3 ln(Za)(m e /m /i ) ln(m M /m e ) and a 2 (Za) 2 ln(Za), respectively. 

The terms (|1^) and ( |17D add up to — 0.496 kHz. The uncertainty due to uncalculated 
terms will be about 0.05 kHz. Including these estimates and using the value of a, i?^ and 



m M /m e from Refs. [18], [|T^| and 



a' 1 = 137.035 997 9 (32) (0.024 ppm) 
Roo = 10 973 731.568 30 (31) m" 1 , 



m,. 



206.768 259 (62), 



(18) 



we find 
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Az/(theory) = 4 463 302.63 (1.34) (0.21) (0.17) kHz , (19) 

where the first and second errors reflect the uncertainties in the measurements of m M and 
oT x listed in (|T8|). The third error is purely theoretical and dominated by the uncertainty 
in the last a(Za) 2 term of (Q). Further reduction of this error is very important and will 
be attempted shortly. The agreement between Az/(theory) and Az/(exp) is excellent, the 
difference being 

- 0.26 (0.16) (1.34) (0.21) (0.17) kHz, (20) 



where the first error is from the experiment (0) and the rest are carried over from fll~9l). 

The result (plj) is obtained using the value of a from (|I~8|) which is determined from the 
quantum Hall effect. Actually, a more accurate value of a is known from the theory and 



measurement of the electron anomalous magnetic moment, which is [fJO 

cT 1 = 137.035 992 22 (94). (21) 
If one uses this instead of (IIH), one finds 



Az/(theory)' = 4 463 303.00 (1.34) (0.06) (0.17) kHz . (22) 

Alternately, if one assumes that the uncertainty due to a is 0.06 kHz and that QED is 
correct to 0.17 kHz, one can determine the muon mass from ([!]) and (^2|). This leads to 

— £ = 206.768 275 (11), (23) 

m e 

which is 5.6 times more accurate than the value quoted in fllTf). This precision is close to 
that expected from the new direct measurement of m M ||. With the new measurement 
of muonium hyperfine structure and further improvement of theory, it will be possible to 
replace ( p3|) by an even better one. Comparison of this result with the directly measured 
m M may be regarded as an alternative way to test the validity of QED. 
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FIGURES 

FIG. 1. Representative diagrams contributing to the a 2 (Za) radiative corrections to the muo- 
nium hyperfine structure in which two virtual photons are exchanged between e~ and The 
muon is represented by x . 

FIG. 2. Two-photon exchange diagrams with fourth-order radiative corrections on the electron 
line. Diagrams which are related to these diagrams by time reversal are not shown explicitly. The 
muon is represented by x . 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402260v2 
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